function [] = fig_sec10_3_backpressureRoutingDual()

close all;

xp = load ('data/backpressureRoutingDual_xp.txt');
qnd = load ('data/backpressureRoutingDual_qnd.txt');
ye = load ('data/backpressureRoutingDual_ye.txt');
objFunc = load ('data/backpressureRoutingDual_objFunc.txt');
jom_xp = load ('data/backpressureRoutingDual_jom_xp.txt');
jom_qnd = load ('data/backpressureRoutingDual_jom_qnd.txt');
jom_ye = load ('data/backpressureRoutingDual_jom_ye.txt');
jom_objFunc = load ('data/backpressureRoutingDual_jom_objFunc.txt');

multiplot ('Time', 'Average number of hops' , objFunc , jom_objFunc , [0 1.3*max(objFunc(:,2))] , '-' , 'Q_{nd}' , qnd , jom_qnd , [0 1.3*max(max(qnd(:,2:end)))] , '-' , 'x_p' , xp , jom_xp , [0 29] , '-');

% multiplot ('Time', 'Network utility' , objFunc , jom_objFunc , [-900 -500] , '-' , 'p_e' , pe , jom_pe , [0 0.2] , '-' , 'u_e' , ue , jom_ue , [0 0.1] , '-');
% multiplot ('Time', 'Network utility' , objFunc , jom_objFunc , [-100 100] , '-' , 'p_e' , pe , jom_pe , [0 1] , '-' , 'u_e' , ue , jom_ue , [0 1] , '-');

jom_objFunc
grad_objFunc = objFunc (end)

max(xp(end,2:end))
max(jom_xp)

max (jom_ye)
max (ye(end,2:end))